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Fundamentals  of  SAR  Polarimetry 

Abstract  A  comprehensive  overview  of  the  basic  principles  of  radar  polarimetry  is  presented.  The  relevant 
fundamental  field  equations  are  first  provided.  The  importance  of  the  propagation  and  scattering  behavior  in 
various  frequency  bands,  the  electrodynamic  foundations  such  as  Maxwell’s  equations,  the  Helmholtz  vector 
wave  equation  and  especially  the  fundamental  laws  of  polarization  will  first  be  introduced:  The  fundamental 
terms  which  represent  the  polarization  state  will  be  introduced,  defined  and  explained.  Main  points  of  view 
are  the  polarization  Ellipse,  the  polarization  ratio,  the  Stokes  Parameter  and  the  Stokes  and  Jones  vector 
formalisms  as  well  as  its  presentation  on  the  Poincare  sphere  and  on  relevant  map  projections.  The  Polarization 
Fork  descriptor  and  the  associated  van  Zyl  polarimetric  power  density  and  Agrawal  polarimetric  phase 
correlation  signatures  will  be  introduced  also  in  order  to  make  understandable  the  polarization  state 
formulations  of  electromagnetic  waves  in  the  frequency  domain.  The  polarization  state  of  electromagnetic 
waves  under  scattering  conditions  i.e.  in  the  radar  case  will  be  described  by  matrix  formalisms.  Each 
scatterer  is  a  polarization  transformer;  under  normal  conditions  the  transformation  from  the  transmitted  wave 
vector  to  the  received  wave  vector  is  linear  and  this  behavior,  principally,  will  be  described  by  a  matrix 
called  scattering  matrix.  This  matrix  contains  all  the  information  about  the  scattering  process  and  the 
scatterer  itself.  The  different  relevant  matrices,  the  respective  terms  like  Jones  Matrix,  S-matrix,  Muller  M- 
matrix,  Kennaugh  K-matrix,  etc.  and  its  interconnections  will  be  defined  and  described  together  with  change 
of  polarization  bases  transformation  operators,  where  upon  the  optimal  (Characteristic)  polarization  states  are 
determined  for  the  coherent  and  partially  coherent  cases,  respectively.  The  lecture  is  concluded  with  a  set  of 
simple  examples. 

4.  Polarimetric  Radar  Optimization  for  the  Coherent  Case 

The  optimization  of  the  scattering  matrices,  derived  for  the  mono-static  case  is  separated  into  two  distinct 
classes.  The  first  one,  dealing  with  the  optimization  of  [5] ,  [G] ,  and  [A] ,  for  the  coherent  case  results  in 
the  formulation  of  ‘ Kennaugh  ’s  target  matrix  characteristic  operator  and  tensorial  polarization  fork  ’  and 
the  associated  renamed  ‘ Huynen  Polarization  Fork’  concept  plus  the  ‘co/cross-polarization  power  density 
plots’  and  the  ‘co/cross-polarization  phase  correlation  plots’,  also  known  as  the  van  Zyl  [79,  71]  and  the 
Agrawal  plots  [78,  90],  respectively,  in  the  open  literature.  The  second  one,  presented  in  Chapter  5,  deals 
with  the  optimization  for  the  partially  polarized  case  in  terms  of  the  ‘ lexicographic  and  the  Pauli-based 
covariance  matrices,  [cj  and  [cj.  respectively’,  as  introduced  in  Sections  3.7  to  3.10,  resulting  in  the 
‘Cloude  target  decomposition  theorems’  and  the  Cloude-Pottier  [27,  57,  58]  supervised  and 
unsupervised  ‘Polarimetric  Entropy  H,  Anisotropy  A,  and  a -Angle  Descriptors’.  In  addition,  the 
‘ polarimetric  contrast  optimization  procedure  ’  dealing  with  the  separation  of  the  desired  polarimetric  radar 
target  versus  the  undesired  radar  clutter  returns  of  which  the  alternate  lexicographic  and  Pauli-based 
covariance  matrix  optimization  procedures  deserve  special  attention  next  to  the  coherent  [.S']  and  partially 
coherent  [A]  matrix  cases. 

4.1  Formulation  of  the  Mono-Static  Radar  Optimization  Procedure  according  to  Kennaugh  for 

the  Coherent  Case 

Kennaugh  was  the  first  to  treat  the  mono-static  polarimetric  radar  optimization  procedure  (see  Fig.  4.1)  for 
optimizing  (3.9)  according  to  the  BSA  formulation 

Ev(r)  =  [5]  E'*  (r)  (4.1) 
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Fig.  4.1  BSA  Optimization  According  to  Kennaugh 


but  with  the  received  field  E'  (r)  being  so  aligned  with  the  incident  field  E'  (r)  with  the  reversal  of  the 
scattered  versus  incident  coordinates  of  the  BSA  system  resulting  in  Kennaugh’s  psuedo-eigenvalue’  [4] 
problem  of 

Opt{[S]}  such  that  [.S’]  E  -  A  E  =  0  (4.2) 


The  rigorous  solution  to  this  set  of  ‘con- similarity  eigenvalue’  problems  was  unknown  to  the  polarimetric 
radar  community  until  the  late  1980’s,  when  Liineburg  [54],  rediscovering  the  mathematical  tools  [116,  117], 
derived  a  rigorous  but  mathematically  rather  involved  method  of  the  associated  con-similarity  eigenvalue 
problem,  not  further  discussed  here,  but  we  refer  to  Liineburg’s  complete  treatment  of  the  subject  matter  in 
[52,  53].  Instead,  here  Chan’s  [77]  ‘ Three-Step  Solution  as  derived  from  Kennaugh’s  original  work  [4],  is 
adopted. 

Three  Step  Procedure  according  to  Chan  [77] 

By  defining  the  polarimetric  radar  brightness  (polarization  efficiency,  polarization  match  factor)  formation 
according  to  (3.26)  and  (3.27)  retaining  the  factor  1/2  (not  contained  in  the  1983  IEEE  Standard  and  in 
Mott’s  textbook)  [76,  102]  as 


2 


[S]  E' 


(4.3) 


in  terms  of  the  terminal  voltage  VR  ,  being  expressed  in  terms  of  the  normalized  transceiver  antenna  height 
h'  and  the  incident  field  E' ,  as  defined  in  Mott  [76]  and  in  [19],  by 


VR=h'T  Es  =hrT[s]E'  with  h  = 
so  that  the  total  energy  density  of  the  scattered  wave  Es  ,  may  be  defined  by 
W  =  Es+Es  =  ([s]E')+([.S]E')=E'+  ([S]+[S])E'  =  E(+[g]E' 

where  [G]  =  [S']1  [S']  defines  the  Graves  power  density  matrix  [7],  first  introduced  by 
Kennaugh  [4,  5]. 


(4.4) 


(4.5) 
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Step  1 

Because  the  solution  to  the  'pseudo-eigenvalue  problem’  of  (4.2)  was  unknown  at  that  time  (1954  until 
1984);  and,  since  [5]  could  be,  in  general,  non-symmetric  and  non-hermitian,  Kennaugh  embarked  instead 
in  determining  the  ‘optimal polarization  states  ’  from  optimizing  the  power  density  matrix  so  that 

(4.6> 

.  .  2  .  .  2 

for  which  real  positive  eigenvalues  v  j  =  \A  ,|  and  v  2  ij  exist  for  all  matrices  [S]  since  [G]  is 
Hermitian  positive  semidefinite  so  that 

v  12  =  ~  | Trace  [a]  ±  \frace2  [g]  -  4Def  [g]]2 1  (4.7) 

where 

v  i  +  v  2  =  invariant  =  Trace{G\  =  Span[S ]  =  [SHH  f  +  \_SHV  f  +  [SVH  f  +  [Svv  f  =  tc4  (4.8) 

v  !  •  v  2  =  invariant  =  Det[G]  =  ( Det[S ]  )(Det[S ]  )*  =  (SHHSVV  -  SHVSVH  ){SHHSVV  -  SHVSVH )'  (4.9) 

For  the  mono-static,  reciprocal  symmetric  [.S] ,  above  equations  reduce  with  SHV  =  SVH  to 

v  i  +  v  2  =  invariant  =  Trace[G]  =  Span[S ]  =  [5^  f  +  l\SHV ]2  +  [Svv  f  =  k 3  (4. 1 0) 

and 

v  j -v  2  =  invariant  =  Det[G]  =  (Det[s])(Det[s]y  =(shhSvv  ~\SHV\2)  (4.11) 

In  order  to  establish  the  connection  between  the  coneigenvalues  of  equation  (4.2)  and  the  eigenvalues  of 
[G]  in  (4.6),  one  may  proceed  to  take  the  complex  conjugate  of  (4.2)  and  insert  back  in  (4.2).  Equation  (4.2) 
has  orthogonal  solutions  if  and  only  if  [.S']  is  symmetric.  The  inverse  step  is  much  more  difficult  to  prove 
and  needs  among  others  the  symmetry  of[<S] ,  which  provides  another  topic  for  future  research. 

As  a  result  of  these  relations,  Kennaugh  defined  the  ‘effective  polarimetric  radar  cross-section  ’  ZK  ,  also 
known  as  ‘ Kennaugh ’s  Polarimetric  Excess  SK^  in  [118],  where 

sK=  Span[S]+  2|/9ef[5']|  (4.12) 

which  comes  automatically  into  play  (also  in  the  present  formulation)  when  representations  on  the  Poincare 
sphere  are  considered,  which  reduces  to  eK  for  the  mono-static  reciprocal  case.  It  plays  an  essential  role  in 

Czyz’s  alternate  formulation  of  the  ‘theory  of  radar  polarimetry’  [110],  derived  from  a  spinorial 
transformation  concept  on  the  ‘generalized polarization  sphere  ’,  being  studied  in  more  depth  by  Bebbington 
[32]. 

Step  2 


RTO-EN-SET-081 


2-3 


Basics  of  SAR  Polarimetry  li 


ORGANIZATION 


Using  the  resulting  solutions  for  vl2  for  the  known  [g]  =  [.S']  [.S’]  and  [5],  the  optimal  transmit 
polarization  states  Ej,pt  and  optimal  scattered  waves  E('ip[|  can  be  determined  as 


^°pti2 


(4.13) 


Step  3 

The  received  optimal  antenna  height  h'opt  is  then  derived  from  (4.4) 


as  h  opt 


(4.14) 


which  defines  the  ‘polarization  match  ’  for  obtaining  maximum  power  in  terms  of  the  polarimetric  brightness 
function  (4.4)  introduced  by  Kennaugh  in  order  to  solve  the  polarimetric  radar  problem  [4]. 

There  exist  several  alternate  methods  of  determining  the  optimal  polarization  states  either  by 
implementing  the  ‘generalized  complex  polarization  p  ’  transformation,  first  pursued  by  Boemer  et  al.  [13]; 
the  ‘con-similarity  transformation  method’  of  Liineburg  [52,  53]  ,  the  ‘spinorial  polarization  sphere 
transformations  ’  of  Bebbington  [32],  and  more  recently  the  ‘Abelian  group  method  ’  of  Yang  [104  -  105].  It 
would  be  worthwhile  to  scrutinize  the  various  approaches,  which  should  be  a  topic  for  future  research. 

4.2  The  Generalized  p  -  Transformation  for  the  Determination  of  the  Optimal  Polarization  States 

by  Boerner  using  the  Critical  Point  Method 

Kennaugh  further  pioneered  the  ‘polarimetric  radar  optimization  procedures’  by  transforming  the 
optimization  results  on  to  the  polarization  sphere,  and  by  introducing  the  co-polarized  versus  cross-polarized 
channel  decomposition  approach  [4]  which  were  implemented  but  not  further  pursued  by  Huynen  [9]. 
Boemer  et  al.  [31,  82],  instead,  proposed  to  implement  the  complex  polarization  ratio  p  transformation  in 
order  to  determine  the  pairs  of  maximum/minimum  back-scattered  powers  in  the  co/cross-polarization 
channels  and  optimal  polarization  phase  instabilities  (cross-polar  saddle  extrema)  by  using  the  ‘critical point 
method ’  pioneered  in  [82].  Assuming  that  the  scattering  matrix  \S(I  IV)]  is  transformed  to  any  other  ortho¬ 
normal  basis  [AB]  such  that 


sAA  SAB 

e'  e' 

~BA  ^  BB 


(4.15) 


with  [G]  given  by  (2.23);  and  SHV  =  SVH  ,  SAB  =  SBA  for  the  mono-static  case,  the  polarimetric  radar 
brightness  equation  becomes 


P  = 


]_ 

2 


(4.16) 


where  the  prime  ‘  refers  to  any  new  basis  [A  B}  according  to  (4.15). 

By  implementation  of  the  Takagi  theorem  [116],  the  scattering  matrix  [.S’  (AB)]  can  be  diagonalized  [52]  so 
that 
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[•s'fis)] 


}  AA 

0 


’  BB . 


A,  ,  0 

0  A  i 


(4.17) 


^  1=  s'aa  (a)  =  (1  +  P  i/7  i)  1  +  2/?  +  P  X2)  exp(2jyi)  =  |2  J  exp(2y'0  ,)  (4.18) 

X  ^(Pi)=  (l  +  AA*)  '(pI2^-  2A  +  *V)  exp(2yy  4)  =  |A2|  exp(2y^  2) 


as  shown  in  [31]. 


(4.19) 


Determination  of  the  Kennaugh  target  matrix  characteristic  polarization  states:  The  expression  for  the  power 
returned  to  the  co-pol  and  cross-pol  channels  of  the  receiver  are  determined  from  the  bilinear  form  to 
become: 

(i)  Cross-pol  Channel  Minima  or  Nulls  (n),  Maxima  (m)  and  Saddle-Optima  (s): 

For  the  Cross-pol  channel  power  Px  with  E'  =  E!  ,  expressed  in  terms  of  the  antenna  length  h 


|2  1  \.  '  J - -  >  |2 


P,  =  \K\  =  2  h  [5]hl  = 


1  { \  *  |2  ’  ' *  o  o  *  '*2  '2  I  o  |2  1  \ 

7  ’  rry  K1 1  \  P  P  ~  ^  iP  2 P  ^  2 \  P  P  ) 

(i+pp-) 


(4.20) 


so  that  for  the  cross-pol  nulls  ( p  xnl  2  ),  for  the  cross-pol  maxima  ( p  xml  2  ),  and  for  the  cross-pol  saddle 
optima  (p  xsll  ),  according  to  the  critical  point  method  introduced  in  [82], 

(4.21) 


^2  n  *  ^ 

4 

^2  2  *  ^ 

P  xnl, 2  O’00  P  xml, 2  —  —  j 

1  2 
Uxj 

=  ±  exp(y(2v  +  n  /2))  p  xsl>2  =  ± 

1  2 

=  ±  exp(y'(2v)) 


where 


i  '  '*  ■%  '  ■% 

P  xn\P  xnl  P  xm\P  xml  P  xslP  xsl 


(4.22) 


which  states  that  there  exist  three  pairs  of  orthogonal  polarization  states,  the  cross-pol  minima  ( p  2  ),  the 
cross-pol  maxima  ( p  xml  2 ),  and  the  cross-pol  saddle  optima  (  p  xsl  2  ),  which  are  located  pair-wise  at 

antipodal  points  on  the  polarization  sphere  so  that  the  lines  joining  the  orthogonal  polarization  states  are  at 
right  angles  to  each  other  on  the  polarization  sphere[82]. 

(ii) Co-pol  Channel  Maxima  ( p  cmX  2 )  and  Minima  or  Nulls  <  p  cnl  2}j_ 

For  the  function  of  the  power  Pc  with  E'  =  E' ,  return  to  the  co-pol  channel  (c) 

111-  1  I  ,7-  ,|2  1 

K  =  tK  =  r  h  r  [S„]  h  =  - — 5 

2  2  (l  +  PP  ) 


(|2  ,  j  +  A  ,1  \p  '2+  A  \A  2p  '2+\A  ,[  p  p 1,2 1 

(4.23) 


the  critical  points  are  determined  from 

P  cmi=P  xni=Q  p'cmi=p'x„i=co »  where  p'cml  p'cm2=~l 


(4.24) 


and  the  co-pol  maxima  are  identical  to  the  cross-pol  nulls  as  was  first  established  by  Kennaugh 
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[4,  5],  and  utilized  by  Huynen  [9].  In  addition  the  critical  points  for  the  co-pol-null  minima  or  nulls  ( p  cn]  2 ) 
are  determined  from  (4.23)  to  be 


P  cnl,2 


rpi.n 

=  + 

2  ) 

M  2IJ 

exp(y'(2v  +  n  / 2)) 


(4.25) 


and  it  can  be  shown  from  above  derivations  that  the  co-pol-null  minima  p  cnX  and  p  cn2 

lie  in  a  plane  spanned  by  the  co-pol-maxima  (cross-pol-minima)  and  the  cross-pol-maxima  and  the  angle 
between  the  origin  of  the  polarization  sphere  and  the  two  co-pol-nulls  is  bisected  by  the  line  joining  the 
orthogonal  pair  of  co-pol-maxima  (cross-pol-minima)  defining  the  target  matrix  critical  angle  2x2y  as  shown 
first  by  Kertnaugh  [4]  leading  to  his  tensorial  polarization  fork  formulation. 

(iii)Orthogonalitv  Conditions  with  Corresponding  Power  Returns: 

The  three  pairs  of  cross-pol-extrema,  the  cross-pol  nulls  ( p  xnl  2  )  being  identical  to  co-pol  maxima 

( p  cm,  2  ),  the  cross-pol  maxima  ( p  xml  2 )  and  cross-pol  saddle  optima  ( p  xsl  2  ),  satisfy  the  orthogonality 

conditions  of  (4.22)  and  (4.24)  which  implies  that  they  are  located  each  at  anti-podal  locations  on  the 
polarization  sphere.  We  note  that  the  co-pol  maxima’  consist  of  one  absolute  maximum  and  an  orthogonal 
local  maximum.  The  corresponding  co/cross-polar  power  returns  become 

Min \P  }  =  PM (p  „, ) =  P„, lp  ,,,  )=  (I  : 

Max{Pr }:  Prw,(p'cml)= l-ij2  P,mXp\„i)=\2  2f  ; 

Max{P,}:  l>Jp„u)=i<|*,f+M2); 

Min{Pc}:  P.(p„,,2)=0; 

(4.26) 

Sad{Px\.  r,(p',„,2)=I(|X,Ml2f) 

The  resulting  co/cross-polar  extrema  are  plotted  on  the  polarization  sphere  shown  in  Fig.  4.2 
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4.3  The  Kennaugh  Target  Characteristic  Scattering  Matrix  Operator,  and  the  Polarization  Fork 

according  to  Boerner 

Kennaugh  was  the  first  to  recognize  that  the  orthogonal  pairs  (Xb  X2)  of  the  cross-pol  nulls  {p  xni2)  or  co' 

pol  maxima  ( pcmX  =  pxnX ,  pcm2  -  pxn2  )  and  the  pair  (  Si,  S2  )  of  cross-pol  maxima  ( p  xml2  )  lie  in  one  main 
cross-sectional  plane  of  the  polarization  sphere  also  containing  the  pair  (  Ci,  C2)  of  non-orthogonal  co-pol 
nulls  ( p cn\  2  ),  where  the  angle  4 \y  between  the  two  co-pol  null  vectors  on  the  Poincare  sphere  is  bisected  by 

the  line  joining  the  two  co-pol  maxima  (cross-pol  nulls).  These  properties  were  first  recognized  explicitly 
and  utilized  by  Kennaugh  for  defining  his  “Spinorial  Polarization  Fork”,  used  later  on  by  Fluynen  to  deduce 
his  ‘Polarization  Fork '  concept. 

Flowever,  Boerner  et  al.  [13,  25,  81,  31,  82],  by  implementing  the  complex  polarization  ratio 
transformation,  were  able  to  relate  the  polarization  state  coordinates  P(p  )  on  the  Polarization  sphere 
directly  to  the  corresponding  p  on  the  complex  polarization  ratio  plane.  Then  according  to  [82],  each 
point  p  of  the  complex  plane  can  be  connected  to  the  ‘zenith  (LC)  ’  of  the  polarization  sphere,  resting 
tangent  to  the  complex  plane  in  its  ‘origin  0  ’  of  the  ‘nadir  (RC)  ’,  by  a  straight  line  that  intersects  the  sphere 
at  one  arbitrary  point ,  where  the  ‘nadir  (RC)  ’  corresponds  to  the  ‘origin  0  ’  of  the  p  -plane,  the  ‘zenith  (z)  ’ 
to  the  p  -circle  at  ‘infinity  (cq  ’  and  the  equator  representing  linear  polarization  states.  Any  two  orthogonal 
polarization  states  are  antipodal  on  the  sphere,  like  ‘zenith  (left-circular)’,  and  ‘nadir  (right-circular)’. 
Utilizing  this  property,  Boerner  and  Xi  [31]  were  able  to  associate  uniquely  three  pairs  of  orthogonal 
polarization  states  at  right  angle  (bi-orthogonal)  on  the  polarization  sphere;  i.e.  the  anti-podal  points  Si,  S2 

( p  xm\  2  )  and  Ti,  T2  ( p  xsl  2  ),  with  SXS2  and  TXT2  being  perpendicular  to  one  another  (bi-orthogonal);  and 
similarly  to  the  line  XxX2  joining  X!  (nadir:  p  xnX=  p  cmX=  0)  and  X2  (zenith:  p  xn2=  p  cm2=  oo  );  where 
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the  co-pol  nulls  ( p  cnX  2  )  lie  on  the  same  main  circle  on  the  complex  plane  of  p  cnl  2  and  p  xnX  2  so  that 

their  corresponding  points  Ci  and  C2  are  symmetric  about  the  diameter  XxX2  which  bisects  the  angle 
between  Ci,  0,  C2  known  as  the  Kennaugh  target  matrix  characteristic  angle  (  2x2  y ).  The  ‘great  cross- 

sectional  plane  ’  containing  SXS2  -L  XxX2  and  Ci,  C2  is  denoted  as  the  ‘Kennaugh  target  matrix 
characteristic  plane’  with  corresponding  great  circle  being  the  ‘Kennaugh  target  matrix  characteristic 
circle’.  Fig.  4.3  define  the  representations  of  the  Poincare  sphere  for  the  general  and  standardized 
polarization  fork  (Fluynen),  respectively,  including  the  proper  definitions  of  ‘Huynen ’s  Geometrical 
Parameters  (cp-target  orientation  angle;  y-target  ship  angle;  z-target  ellipticity  angle,  p  —  tan  a  extp(jS) 
next  to  ‘Kennaugh  ’s  target  matrix  characteristic  angle  y  ’. 


This  concludes  the  description  of  the  Kennaugh  polarimetric  target  matrix  characteristic  operator’-,  which 
was  coined  ‘the  polarization  fork’  by  Fluynen  [9]. 

4.4  Huynen’s  Target  Characteristic  Operator  and  the  Huynen  Polarization  Fork 

Fluynen  [9],  utilizing  Kennaugh’s  prior  studies  [4,  5],  elaborated  on  polarimetric  radar  phenomenologies 
extensively,  and  his  “Dissertation  of  1970:  Phenomenological  Theory  of  Radar  Targets”  [9],  re-sparked 
international  research  on  Radar  Polarimetr y,  commencing  with  the  studies  by  Poelman  [10,  11],  Russian 
studies  by  Kanareykin  [122],  Potekhin  [123],  and  others  [1]. 

Fluynen  cleverly  reformulated  the  definition  of  the  polarization  vector,  as  stated  in  [9],  so  that  group- 
theoretic  Pauli-spin  matrix  concepts  may  favorably  be  applied  which  also  serve  for  demonstrating  the 
orientation  angle  invariance  which  Fluynen  coined  ‘de-psi-ing  (de-y/  -ingf  using  1//  for  denoting  the 
relative  polarization  ellipse  orientation  angle.  Flere,  we  prefer  to  divert  from  our  notation  by  rewriting  the 
parametric  definition  of  the  polarization  vector 


P(l  E  \,(j),y/,x)  =|  E  |  exp(y^) 


cos^ 

sin^z 


-sin^z 

cosyz 


cos  ^ 
j  sin  X 


as 


p(a,  a,  (f>,  r)  =  a  exp {ja) 


costp 

1 — 

(/i 

1 

cosr 

’c n 
_ 1 

cosf 

-  j  sin  r 

(4.27) 


which  with  the  use  of  the  Pauli-spin  matrices  [cr  ]  defined  in  (2.14),  the  Fluynen  quaternion  group 
definitions  may  be  re-expressed  as  [  /  ]  =  [  cr0  ] ,  [ »/  ]  =  — y  [  cr3  ] ,  [if]  =  y[cr2],  [  Z.  ]  =  — y  [  CTj  ] 


p  {a,a,(j),r)  =  a  exp  {ja)  exp(^[  J])  exp(  r[K] ) 


1 

0 


(4.28) 


with  exp(^[J])  =  cos^[/]  +  sin^[y]  and  exp(r[K])  =  cos  r[/]  +  sin  r[K] 


In  this  notation  the  orthogonal  polarization  vector  p  becomes 


p±(a± ,a±,<f>,r)  =  a±  exp(y a±) exp {(f  +  — )[ J] } exp(-r[K]) 


(4.29) 
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so  that  pp*=a2,  pp^=0. 

Utilizing  this  notation,  the  transformed  matrix  [.S’  (AB)\  becomes 

[S'(AB)]  =  [U]t[S(HV)][U]  (4.30) 


with  the  orthonormal  transformation  matrix  [17]  defined  in  (2.23),  which  may  be  recasted  with  m  =  pm 
the  so-called  maximum  or  null  polarization  as  defined  in  [31],  into 

[t/]  =  [mmj  (4.31) 

Because  of  the  orthonormal  properties  of  m  and  m  ,  which  satisfy  the  con-similarity  eigenvalue  equation 
[82],  the  off-diagonal  elements  of  [S  ( AB )]  vanish.  This  in  turn  can  be  used  to  solve  for  p  in  (2.23),  and 
hence  for  m  and  m  ,  without  solving  the  consimilarity  eigenvalue  problem  of  (4.6).  The  complex 
eigenvalues  pxnl  2 ,  defined  in  (4.26)  are  renamed  as  5, 2  and  were  defined  by  Huynen  as 


x,  =  m  exp  {2  j(u  +  /?)}  s2  =  m  tan2  y  exp  {-2  j{u  -ft)}  (4.32) 

so  that  [5  (. AB )]  of  (  4.30)  becomes 


[SW)]  =  [£r(m,m±)] 


mexp{2j(u  +  /3)}  0 

0  m  tan2  y  exp  {-2  /( v  -  /?)} 


[(7*(m,mx)]7 


(4.33) 


where  m  =  crK,y,i//,Tm,u,  and  /3  are  the  Huynen  parameters,  and  in  denoting  the  target  matrix 
magnitude,  may  be  identified  to  be  “Kennaugh’s  polarimetric  excess <JK”  defined  in  (  4.12);  and 
m(i//,  Tm  )  may  be  re-normalized  as 


cos  y/ 

-sin  y/ 

COS  X m 

=  exp(i//[j])exp(rm[K]|^ 

sin^z 

cos  yr 

_-7'sinr)H 

(4.34) 


and  finally  with 


mcxp{2j(u  +  /])) 
0 


m 


tan2  y  cxp{- 2j(v - /()} 


=  m  exp(//?) 

"1  0 

r\  .  2 

0  tan  y 

ex 


pHl]) 


(4.35) 


and  [Z]  representing  the  third  modified  Pauli-spin  matrix,  satisfying  in  Huynen’ s  notation 


[lHj]Ik]  =  -Ik][j]- 


-j  0 
0  j 


W2=-M 


(4.36) 


we  obtain  Huynen’ s  target  matrix  characteristic  operator 


[H(m,  y,  p;  (j)m ,  u,  xm )]  =  [t/*(i|/,  xm ,  o)] 


m 


1 


0 


0  tan1  y 


[f/*(i|/,Tm,u)]/exp(yP) 


(4-37) 
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where 

rm ,  u)]cxp(i//[j  ])exp(r„,  [k])cxP(u  M)  (4.38) 

representing  the  Eulerian  rotations  with  2 y/,  2rm,  2u  about  the  bi-orthogonal  polarization  axes 
S{S2  (connecting  the  two  cross-pol  maxima),  XxX2  (connecting  the  two  cross-pol  nulls,  or  equivalently,  co- 
pol  maxima),  and  TXT2  (connecting  the  two  saddle  optima),  respectively,  with  more  detail  given  in  Boemer 
and  Xi  [31].  Huynen  pointed  out  the  significance  of  the  relative  target  matrix  orientation  angle  (P  =  (j>  —  \f>  , 
where  (f)  denotes  the  antenna  orientation  angle,  and  that  from  definition  of  m(i//,  rm  )  in  (4.34),  it  can  be 

shown  that  it  can  be  eliminated  from  the  scattering  matrix  parameters  and  incorporated  into  the  antenna 
polarization  vectors  (‘ de-psi-ing :  de-y/  -ing’)  ,  and  that  the  Huynen  parameters  are  orientation  independent, 
which  was  more  recently  analyzed  in  depth  by  Pottier  [58].  The  Eulerian  angle  are  indicators  of  a  scattering 
matrix’s  characteristic  structure  with  o  denoting  the  so-called  ‘ skip-angle ’  related  to  multi-bounce 
scattering  (single  versus  double),  Tm  denotes  the  helicity-angle  and  is  an  indicator  of  target  symmetry 
Tm  =  0  or  non-symmetry,  and  J3  is  the  absolute  phase  which  is  of  particular  relevance  in  polarimetric  radar 
interferometry. 


[S]  =  [<r(p)]exp(v[I]*)m 


Mp)Y 


1  0 
0  tan2y 


exp(v[L]*)r[CT(p)]rexp(/§) 


,-J** 


«  ■  I  *i  I 

/  „  /  /  r/  / 

P*/i1  =  "  ”  PcjbI  >  Pxn2  "  '  “  Pcn2 

Pf»i,2  “  *  tan(*/2  -  y  )exp[/(2y  +  ic/2)] 

Pimi.2 =  ±exp[/'(2Y  +«/2)]  (LR) 

P;j1i2  =  ±expO'2Y)  (45° /1 35° ) 

Huynen's  Solution  of  the  Polarization  Fork 

'1  0 


i , — /  ft 


[y*;(’l'.tw,v)]expOO 

0  tan2Y. 

Fig.  4.3  Huynen’s  Polarization  Fork  (Xi-Boemer  Solution) 


4.5  Alternate  Coherent  Scattering  Matrix  Decompositions  by  Kroggager  and  by  Cameron 

Another  class  of  scattering  matrix  decomposition  theorems  [124,  30]  were  recently  introduced,  and  are  also 
expressed  in  terms  of  the  Pauli  spin  matrix  sets  i  =  0, 1,  2,  3),  by  associating  elementary  scattering 

mechanisms  with  each  of  the  [cr  ] ,  so  that  for  the  general  non-symmetric  case 
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\s{a,b,c,dy\ 


a  +  b 
c  +  jd 


c-jd 

a-b 


=  «[°o  ]  +  bWi  ]  +  cWi  ]  +  d[(Ty  ] 


(4.39) 


where  a,  b,  c  and  d  are  all  complex.  Above  ‘ coherent ’  decomposition  may  be  interpreted  in  terms  of  four 
‘ elementary  deterministic  point  target  ’  scattering  mechanisms,  viewed  under  a  change  of  wave  polarization 
basis,  where 

a  -  corresponds  to  single  scattering  from  a  sphere  or  plane  surface 

b  -  corresponds  to  di-plane  scattering  (4.40) 

c  -  corresponds  to  di-plane  scattering  with  a  relative  orientation  of  45° 

d  -  corresponds  to  anti-symmetric  ‘helix-type’  scattering  mechanisms  that  transform  the  incident  wave 
into  its  orthogonal  circular  polarization  state  (helix  related) 

Krogager  [30]  and  Cameron  [124,  125],  among  others,  in  essence  made  use  of  this  decomposition  for  the 
symmetric  scattering  matrix  case  by  selecting  the  desirable  combinations  of  the  [cr  ]  that  suits  their  specific 
model  cases  best. 

In  the  ‘ Krogager  Approach  ’,  a  symmetric  matrix  [.S),.  ]  is  decomposed  into  three  coherent  components, 
which  display  the  1 physical  meaning’  of ‘sphere’,  ‘diplane’,  and  ‘helical  targets’,  where 

[SK  [a,b,d)\  =  a  [Ssph ]  +  //  exp (/>/;)  [Sdi ]  +  rj  cxp(/7//)  [She, ]  (4.41) 


with  additional  direct  associations  with  the  Pauli  matrices  defined  in  (2.14)  given  by 


1  0 
0  1 


1  0 

0  -1 


0  ~j 

j  o 


(4.42) 


This  decomposition  is  applied  directly  to  the  complex  [5]  matrix  imagery,  and  results  into  a  rather  efficient 
sorting  algorithm  in  terms  of  the  three  characteristic  ‘ feature  sorting  base  scatter  images’.  Using  color 
composite  presentations  for  the  three  classes  then  allows  for  the  associated  ‘unsupennsed  feature  sorting  ’. 
This  feature  sorting  method  has  been  applied  rather  successfully  in  the  interpretation  of  various  geo- 
environmental  (forestry,  agriculture,  fisheries,  natural  habitats,  etc.)  as  well  as  in  law  enforcement  and 
military  applications. 


In  the  ‘ Cameron  Approach’,  the  matrix  \Sc(a,b,c,d)\  is  decomposed,  by  separating  the  non¬ 
reciprocal  [s^  ]  from  the  reciprocal  component^  ]  =  [sc  ]  via  an  orientation  angle  (j)' ,  and  by  further 


decomposing  the  latter  into  two  further  components, 
the  angle  z  ,  so  that 


nmax  j 

oc  and 


^rmin 


\Sc{a,b,c,d^\  =  a  cos^'  jeosr 


+  sin  z 


5;n;;;  |+ a  sin  </>'[sc J 


,  with  linear  eigen-polarizations  via 


(4.43) 


which  is  further  analyzed  in  [40]  and  its  limitations  are  clearly  identified  in  [27]. 

Of  the  many  other  existing  [*S («, Z?,  c,  c/ )]  matrix  decomposition  theorems,  mostly  derived  from  alternate 
formulations  (4.39)  of  the  Pauli  spin  matrix  set  ipp  ([cr  ] ,  i  =  0, 1,  2, 3),  defined  in  (2.14),  the  three  examples 

of  the  Huynen,  the  Krogager  and  the  Cameron  decompositions,  it  becomes  apparent  that  there  exists  an 
infinum  of  decompositions  non  of  which  is  unique  and  all  of  them  are  basis  dependent  and  require  a  priori 
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information  on  the  scatterer  scenario  under  investigation.  Yet  for  specific  distinct  applications  all  of  them 
may  serve  a  useful  purpose  which  highly  superior  to  any  non-polarimetric  or  partially  polarimetric  treatment. 


(a)  San  Francisco  Image  (b)  SDH  Decomposition 


Fig.  4.4  Polarimetric  Decompositions:  Krogager,  1993.  (a)  Original  San  Francisco  POLSAR 
image  with  RGB  color  coded  by  |HH-VV|,  |HV|  and  |HH+VV|,  respectively. 

(b)  Sphere  (Blue),  diplane  (Red)  and  helix  (Green)  decomposition  (SDH)  decomposition. 

4.6  Kennaugh  Matrix  Decomposition  of  Huynen’s  Matrix  Vector  Characteristic  Operator 

The  ‘ Kennaugh  target  matrix  characteristic  operator ’  can  also  be  derived  from  the  Kennaugh  matrix,  as  was 
shown  by  Boemer  et  al.  [78,  90,  31,  82];  using  the  Lagrangian  multiplier  method.  A  more  elegant  method 
was  recently  devised  by  Pottier  in  order  to  highlight  the  importance  of  Huynen’s  findings  on  the  ‘ target 
orientation  i//  invariance  ’  for  both  the  Kennaugh  and  the  Lexicographic  Covariance  matrix  representations; 
and  most  recently  by  Yang  [119  -  121],  who  shed  more  light  into  the  properties  of  the  ‘  equi-power -loci’  as 
well  as  ‘equi-correlation-phase-loci  ’  which  deserve  careful  future  attention  but  will  not  be  further  analyzed 
here. 

Instead,  we  return  to  Huynen  [9],  who  provided  further  phenomenological  insight  into  the  properties  of  the 
Kennaugh  matrix  [if],  by  redefining  [,S]  for  the  symmetric  case  in  terms  of  the  limited  set  of  Pauli  spin 
matrices. 


[.S(a,b,c)]  = 


a  +  b  c 
c  a-b 


=  a[a0\  +  b[al]  +  c[a2] 


(4.44) 


so  that  with  the  formal  relation  of  [if]  with  [.S'] ,  obtained  via  a  Kroenecker  product  multiplication  as 

[k]  =  2[a]t'  (4.45) 

insertion  of  (4.44)  into  (4.45)  yields,  using  Huynen’s  notation 
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(A0+B0) 

Fv 

H r 

( ~A0+B0 ) 

D r 

(A0  +  B{ ) 

E ; 

H r 

D> r 

E* 

{a-b«) 

where 


(4.46) 


\a\2=2A0  ab*  =  C-jD 


\b\~  =  B0+B{  be  =  E  +  jF  (4.47) 

\c\~=B0-Bl  ac*=H  +  jG 

Recognizing  that  H  becomes  zero  for  proper  ‘de-y/  -ing’  of  [.S']  as  defined  in  [//]  of  (4.37)  by  removing 
the  iff  rotational  dependence,  and  inserting  it  into  the  antenna  descriptors,  he  was  able  to  redefine  his 
Kennaugh  matrix  coefficients  such  that 
H¥=Csin2y/  C¥  =Ccos2y/ 

Gv=Gcos2y/-Dsm2y/  Ev  =  E  cosAy/ +  F  smAy/  (4.48) 

D  =  G  sin  2y/  +  D  cos  2y/  F  =  —E  sin  Ay/  +  F  cos  Ay/ 


so  that  the  ‘ de-y /  -ed’  [A]  becomes,  by  removing  the  ^-dependence  from  [A  ]  and  incoiporating  it  into 
the  antenna  polarization  and  Stokes’  vectors,  respectively,  such  that 


(A0  +  Bn ) 

F 

C 

o 

II 

F 

(-  Aq  +  Bn ) 

G 

D 

C 

G 

(A0  +  B1 ) 

E 

o 

ll 

D 

E 

o 

cq 

1 

O 

cs 

and  expressed  in  terms  of  the  Huynen  parameters 


A0=Qf  cos2  2rm 
B0  =  q( 1  +  cos2  2 y-f  cos2  2 rm ) 
C  =  2Q  cos  2  y  cos  2vm 
F  =  2Q  cos  2  y  sin  2rm 
Q  =  m2(8cos4  /)  * 


G  =  Qf  sin4rm 

5,  =g(l  +  cos227-/(l  +  sin22rm)) 
D  =  Q  sin2  2y  sin  Ao  cos  2 zm 
E  =  -0sin2  2y  sin4t>sin2rm 
/  =  1-sin2  2y  sin2  so 


(4.49) 


(4.50) 


Huynen’ s  decomposition  was  greeted  with  comprehension,  steadily  but  slowly  so  are  his  phenomenological 
argumentations.  We  note  that  its  uniqueness  is  not  guaranteed,  because  it  is  not  basis-independent  as  was 
shown  by  Pottier  [126]. 

4.7  Optimization  of  the  Kennaugh  matrix  \K\  for  the  Coherent  Case  Using  Stokes’  Vector 

Formulism 

Using  the  Lagrange  multiplier  method  applied  to  the  received  power  matrices  for  the  mono-static  reciprocal 
case  in  terms  of  the  Kennaugh  (Stokes  reflection)  matrices  [^v]>  and  [Km  ]  of  (3.30)  and  (3.31), 
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respectively,  derived  in  [113],  enables  one  to  determine  the  characteristic  polarization  states  similar  to  the 
generalized  p-transformation  methods. 

For  simplicity,  the  transmitted  wave  incident  on  the  scatterer  is  assumed  to  be  a  completely  polarized  and 
normalized  wave  q  so  that 


and  re-stating  the  received  power  expression  defined  in  (3.26),  (3.29),  and  (3.34)  as  functions  of  the  optimal 
Stokes  parameters,  where 


p,  =  <f  Ki  q'  =  p<  p,  =  <f  A]  q'  =  p,  fit 

p,  =  q"  kk  =p„ 


are  subject  to  the  constraint  of  (4.51). This  requirement  dictates  the  use  of  the  method  of  Lagrangian 
multipliers  to  find  the  extrema  of  the  received  powers  P  ,  Px ,  and  Pm  .  Reformulating 
the  equation  of  constraint  to  be  given  by 

<t>  (q\ ,  q*2,  q\  )=  (q\2  +  q2  +  qff1  -\  =  i)  (4.53) 

then  the  Lagrangian  multipliers  method  for  finding  the  extreme  values  of  any  of  the  three  returned  power 
expressions/)  {c[[ ,  q‘2 ,  q\ )  results  with  /  =  c,m,x  in 


dPe  d(j) 
dq\  * dq\ 


i=l,  2, 3 


(4.54) 


For  the  corresponding  ‘’degenerate  deterministic ’  (purely  coherent)  Kennaugh  matrix  [A'],  insertion  of  the 
corresponding/)  {c/[,  q‘2,  q'^)  into  (4.54)  results  in  a  set  of  Galois  equations  yielding  for: 

(i)  the  extreme  ‘co-polar  channel  power  P  ’  four  solutions:  two  maxima  pcm ,  2  which  are 

orthogonal,  and  two  minima  (nulls)  which  may  in  an  extreme  pathological  case  be 
orthogonal  (or  identical  but  generally  are  not);  and  not  one  more  or  not  one  less  of  any 
of  these  extrema; 

(«)  the  returned  ‘cross-polar  channel  power  Px  ’  with  six  extreme  solutions,  being  the 
three  non-identical  pairs  of  orthogonal  polarization  states:  the  cross-polar  maxima 
Pxm\  2  ’  cross-polar  minima  pxnX  2 ,  and  the  cross-polar  saddle  optima  pxsX  2 . 
and  not  one  more  or  not  one  less  of  any  of  these  extrema; 

(iii)  the  returned  power  for  the  ‘matched  antenna  case  Pm  ’  yields  only  two  solutions 
being  identical  to  the  ‘co-polar  maxima  pair  pcml  2  =  pxnX  2  ’  and  not  one  more  or 
not  one  less  of  any  of  these  extrema. 

In  summary,  Opt  ]Pc  (q'x,  q‘2,  yields  always  exactly  four  solutions;  Opt  \PX  )}  yields  always 

exactly  six  solutions  (or  equivalently  three  “bi-orthogonal”  pairs  of  orthogonal  polarization  states); 
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and  Opt  ]Pm  \Cj[ ,  cp ,  cf  j]  yields  always  two  solutions  (or  equivalently,  exactly  one  pair  of  orthogonal 

polarization  states).  This  represents  an  important  result  which  also  holds  for  the  partially  polarized  case 
subject  to  incidence  of  a  completely  polarized  wave. 

A  comparison  of  methods  for  a  coherent  scatterer  of  a  specifically  given  Sinclair  matrix  [.S'] ,  with 
corresponding  [G]  and  [A!],  plus  \Kc\,  \Kx\,  and  \Km ] ,  analyzed  in  [1 13],  clearly  demonstrates  that  all  of 
the  methods  introduced  are  equivalent. 

4.8  Determination  of  the  Polarization  Density  Plots  by  van  Zyl,  and  the  Polarization  Phase 

Correlation  Plots  by  Agrawal 

In  radar  meteorology,  and  especially  in  ‘Polarimetric  Doppler  Radar  Meteorology  ’  the  Kennaugh  target 
matrix  characteristic  operator  concept  was  well  received  and  further  developed  in  the  thesis  of  Agrawal  [78]; 
and  especially  analyzed  in  depth  by  McCormick  [127],  and  Antar  [128]  because  various  hydro-meteoric 
parameters  can  directly  be  associated  with  the  Huynen  or  alternate  McCormick  parameters.  In  radar 
meteorology,  the  Poincare  sphere  visualization  of  the  characteristic  polarization  states  has  become 
commonplace;  whereas  in  wide  area  SAR  remote  sensing  the  co/cross-polarization  and  Stokes  parameter 
power  density  plots  on  the  unwrapped  planar  transformation  of  the  polarization  sphere  surface,  such  as 
introduced  independently  -  at  the  same  time  -  in  the  dissertations  of  van  Zyl  [79]  and  Agrawal  [78],  are 
preferred. 


Because  of  the  frequent  use  of  the  ‘co/cross-polarization  power  density  plots’,  P(p),Pi  (p)  and  Px(p); 
and  the  equally  important  but  hitherto  rarely  implemented  ‘co/cross-polarization  phase  correlation  plots 
Rc(p),  Rcl(p)  and  Rx(p),R  (p)  ’,  those  are  here  introduced.  Following  Agrawal  [78],  who  first 
established  the  relation  between  the  ‘  Scattering  Matrix  Characteristic  Operators  of  Kennaugh  and  Huynen ’ 
with  the  ‘polarimetric  power-density/ phase-correlation  plots  ’,  we  obtain  for  the  reciprocal  case  SAB  =  SBA 


[QJ  =  <  f; 


3  L  A3  L 
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sa 


Ms  A 

-t  \2 


'  BB 


(4.55) 


re-expressed  in  terms  of  the  co/cross-polarimetric  power  density  expressions: 


PXp)  =  {\Saa\1)  P,l(p)  =  {\Sbb  I2)  P,(P)  =  {  (4.56) 


and  the  co/cross-polarization  phase  correlation  expressions: 

K(p)=(s.A)  PcXp)=(sbbs'aa)  RAp)={sA)  KAp)={s,A)  (4.57) 


so  that  [Ci3(/?)]  may  be  rewritten  according  to  (3.68)  -  (3.72)  as 


[c!I  Ml 


p.(p) 
4i  r,(p)‘ 

rM 


H  R,(p)  R,(p) 

2P,(p)  RrUp)’ 

PrHp)  R(p) 


(4.58) 


satisfying  according  to  (3.70)  and  (3.71)  the  following  inter-channel  and  symmetry  relations 
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P±  = 
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--1/ 


P  ) 


=  Pc  (P) 


=  P,{P) 
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P  ) 


=I^W 


(4.59) 


Pi  = 


--1/ 


P  J 


=  |*»| 


and  frequently  also  the  degree  of  polarization  Dp(p)  and  the  degree  of  coherency  //( p)  in  terms  of  the 
directly  measurable  Pc(p),PcL(p)  and  PJp) ;  Rc(p),  Rc±(p)  and  Rx(p),  Rx±(p)  ,  provided  a 

‘dual-orthogonal,  dual-channel  measurement  system  for  coherent  and  partially  coherent  scattering  ensembles 
is  available  requiring  high-resolution,  high  channel  isolation,  high  side-lobe  reduction,  and  high 

sensitivity  polarimetric  amplitude  and  phase  correlation,  where 


Mp)= 


■JpJPPp} 


>  o.(p)= 


k(p)-pp)f  +4K(p)| 

(Pjp)+ P„(p)) 


2  1/2 


,  where  0  <  p(p)  <  Dp  (p)  <  1 

(4.60) 


and  for  coherent  (deterministic)  scatterers  p  —  Dp  =  1 ,  whereas  for  completely  depolarized  scatterers 
P  =  Dp  =  0. 


The  respective  power-density  profiles  and  phase-correlation  plots  are  then  obtained  from  the  normalized 

/  \  J  L  J  L  J L  TT 

polarimetric  radar  brightness  functions  as  functions  of  ( (/),  r)  with  — <(/)<—,  — —  <  r  <  —  so  that 

Vaa(0>  r)  =  pr(0’r)[S(RTV)]p(0,r) 

Pab(0>  t)  =  pf±(0>  T)[S(HV)]p(<f>,  t)  (4.61) 

y ba  (M  =  p  r  (</>,  r)[S(//F)]p±  (<f>,  r) 

V m  (M  =  P 1  (M  [S  (HV)]  p  x  (</>,  r) 

where 

Pc  =1  V AA  |2=l  $Aa(0’  r)  I' 

Px=\Vab\2=\Sab^)\2 

Rc  =1  <I>AA  ~  <t>BB  1=1  arg  VAA  (^  T)  ~  afg  V SB  (0>  T)  I  (4‘62) 

R,  (0,  t)  =|  <t>AA  -  (pAB  |=|  arg  VAA  (</>,r)~  arg  V AB  (</>,t)\ 

R,1  T)  =1  hs  ~  BA  1=1  arg  y SB  (</>’*)-  arg  y BA  r)  I 


In  addition,  the  Maximum  Stokes  Vector  <\()MAX  ,  and  the  maximum  received  power  density  Pm  may  be 
obtained  from 


P,„ (<f>,  r)  =  q0(</>,  r)  =  [K]q{(j),  r) 


(4.63) 


where  examples  are  provided  in  Figs.  4.5  and  4.6  for  one  specific  matrix  case  [31,  82]  given  by 


[S(tfT)]  = 


V 
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~j 


(4.64) 
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(a)  The  Kennaugh  spinorial  (Huynen) 
polarization  fork 


LC 


(b)  Associated  optimal  polarization  states 


(c)  Power  density  Co-pol 


(d)  Power  density  x-pol 


R#l»bv*  Co-poianMkon  phaaa  plot 


(e)  Phase  correlation  co-pol  (f)  Phase  Correlation  x-pol. 


Fig.  4.5  Power  Density  and  Phase  Correlation  Plots  for  eq.  4.64  (by  courtesy  of  James  Morris) 

4.9  Optimal  Polarization  States  and  its  Correspondence  to  the  Density  Plots  for  the  Partially 

Polarized  Cases 

According  to  the  wave  dichotomy  portrayed  for  partially  polarized  waves,  there  exists  one  case  for  which  the 
coherency  matrix  for  the  partially  polarized  case  may  be  separated  into  one  fully  polarized  and  one 
completely  depolarized  component  vector  according  to  Chandrasekhar  [34].  This  principle  will  here  be 
loosely  applied  to  the  case  for  which  a  completely  polarized  wave  is  incident  on  either  a  temporally 
incoherent  (e.g.,  hydro-meteoric  scatter)  or  spatially  incoherent  (e.g.,  rough  surface  viewed  from  different 
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depression  angles  as  in  synthetic  aperture  radar  imaging).  This  allows  us  to  obtain  a  first  order  approximation 
for  dealing  with  partially  coherent  and/or  partially  polarized  waves  when  the  polarimetric  entropy  is  low;  for 
which  we  then  obtain  the  following  optimization  criteria: 
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Fig.  4.6  Power  Density  and  Phase  Correlation  Plots  for  eq.  4.64  (by  courtesy  of  James  Morris) 


The  energy  density  arriving  at  the  receiver  back-scattered  from  a  distant  scatterer  ensemble  subject  to  a 
completely  polarized  incident  wave  may  be  separated  into  four  distinct  categories  where  the  Stokes  vector  is 
here  redefined  with  q  ;)  and  q;(  denoting  the  completely  polarized  and  the  unpolarized  components, 
respectively 


DPq  o 

_  1 
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ST 

o 

_ 1 

<h 
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= 

+ 

0 

q2 

q2 

*7 3 

0 

(4.65) 


and  q  as  well  as  D  were  earlier  defined,  so  that  the  following  four  categories  for  optimization  of  partially 
polarized  waves  can  be  defined  as 


q0  total  energy  density  in  the  scattered  wave  before  it  reaches  the  receiver 

q0Dp  completely  polarized  part  of  the  intensity 
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q0  (1  -  Dp  )  noise  of  the  unpolarized  part 

(4.66) 

—  q0  (1  +  Dp  )  maximum  of  the  total  receptable  intensity,  the  sum  of  the  matched  polarized 

Considering  a  time-dependent  scatterer  which  is  illuminated  by  a  monochromatic  (completely  polarized 
wave)  E' ,  for  which  the  reflected  wave  E '  is,  in  general,  non-monochromatic;  and  therefore,  partially 
polarized.  Consequently,  the  Stokes  vector  and  Karnaugh  matrix  formulism  will  be  applied  to  the  four  types 
of  energy  density  terms  defined  above  in  (4.66). 

4.10  Optimization  of  the  Adjustable  Intensity  D pq0 

The  energy  density Dpq0,  contained  in  the  completely  polarized  part  q  of  q,  is  called  the  adjustable 

intensity  because  one  may  adjust  the  polarization  state  of  the  receiver  to  ensure  the  polarization  match  as 
shown  previously  for  the  coherent  case.  We  may  rewrite  the  scattering  process  in  index  notation  as 

3 

q':  -  X  K  ,qti  where  j  =  0,1, 2, 3  (4.67) 

j= o 


part  and  one  half  of  the  unpolarized  part:  \d  pq{]  ]+  I  —  (i  — 


q  o 


The  adjustable  intensity  D  q0  can  be  re-expressed  as 
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(4.68) 


where  the  q\  are  the  elements  of  the  Stokes  vector  of  the  transmitted  wave.  The  partial  derivative  of 
(Dpqoy  with  respect  to  q[  can  be  derived  as 


d(Dpqs0)2 

dq‘k 


3  3 


i  =  1  vqk  i=l  i= 1  j= 0 


(4.69) 


For  optimizing  the  adjustable  intensity,  we  apply  the  method  of  Lagrangian  multipliers,  which  yields 


8(Dp^J 

K  dq‘k 


i= 1  y=0 


where  (/)  is  the  constraint  equation 

Aq\ » ^ » ^ ) = fc2  ’  v12  >  Vi2  Y2  - 1 = 0 


(4.70) 


(4.71) 


Equation  (4.70)  subject  to  (4.71)  constitutes  a  set  of  inhomogeneous  linear  equations  in 
q\  C4  and  q‘i(ju),  with  solutions  as  three  functions  of/u.  Substituting  q‘ {ju;i  =  1, 2, 3)  into  the 

constrained  condition  (4.71)  leads  to  a  sixth-order  polynomial  Galois  equation  of//.  For  each  p  value, 
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q[,  q2,  q\,  and  Dpqs0  are  calculated  according  to  (4.67)  to  (4.69).  The  largest  (or  smallest)  intensity  is  the 
optimal  intensity;  the  corresponding  q  is  the  optimal  polarization  state  of  the  transmitted  wave. 


4.11  Minimizing  the  Noise-Like  Energy  Density  Term:  qs  (l  —  D  ) 

An  unpolarized  wave  can  always  be  represented  by  an  incoherent  sum  of  any  two  orthogonal  completely 
polarized  waves  of  equal  intensity  [14,  15], which  leads  to  50%  efficiency  for  the  reception  of  the  unpolarized 
wave.  In  order  to  receive  as  much  ‘polarized  energy’  as  possible,  the  noise-like  energy  needs  to  be 
minimized.  The  total  energy  density  of  the  unpolarized  part  of  the  scattered  wave  is  given  by: 


(l--D,k=  <-Dr9i=tK<,i9', 


7=0 


1 


I ITA 

i= 1  V7=0 


(4.72) 


Hitherto,  no  simple  method  was  found  for  finding  the  analytic  closed  form  solution  for  the  minimum; 
instead,  numerical  solutions  have  been  developed  and  are  in  use. 


4.12  Maximizing  the  Receivable  Intensity  in  the  Scattered  Wave:  —  q0(l  +  Dp) 


The  total  receivable  energy  density  consists  of  two  component  parts:  100%  reception  efficiency  for  the 
completely  polarized  part  of  the  scattered  wave  and  50%  reception  efficiency  for  the  unpolarized  part.  The 
resulting  expression  for  the  total  receivable  intensity: 


T  (l  +  D„  k  =  Drq-„  +  \  (l  -  B,  k  =  \  £  K.fl)  +  \ 


7=0 


2V 


3  3 


V  «=1  7=0 


(4.73) 


can  only  be  solved  using  numerical  analysis  and  computation.  The  resulting  maximally  received  Stokes 
vector  is  plotted  in  Fig.  4.7  (where  q  was  replaced  by  p );  and  we  observe  that  for  the  fully  polarized  case  no 
‘ depolarization  pedestal ’  exists.  It  appears  as  soon  as  p  <  1 ,  and  for  p  =  0  it  reaches  its  maximum  of  0.5  for 
which  the  polarization  diversity  profile  has  deteriorated  into  the  ‘ flat  equal  power  density  profile",  stating 
that  the  1 polarization  diversity ’  becomes  meaningless. 


p=l 
coherent 
point  scatterer 


p  =  .8 

distributed  partially 
coherent  scatterer 


p  =  0 

total  polarization 
noise 


Fig.  4.7  Optimal  Polarization  States  for  the  Partially  Polarized  Case 
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In  conclusion,  we  refer  to  Boemer  et  al.  [31,  82],  where  an  optimization  procedure  for  (4.65  -  4.68)  in  terms 
of  [/<"]  for  a  completely  polarized  incident  wave  is  presented  together  with  numerical  examples.  It  should  be 
noticed  here  that  Yang  more  recently  provided  another  more  elegant  method  in  [1 19  -  121]  for  analyzing  the 
statistical  optimization  procedure  of  the  Kennaugh  matrix. 
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6.  Appendices 

A.  The  Standard  Kronecker  Tensorial  Matrix  Product 

Consider  a  matrix  [A]  =  [#y  ]  of  order  ( mxn )  and  a  matrix  [5]  =  J  of  order  (rxs).  The  Kronecker  product 

of  the  two  matrices,  denoted  [A]  ®  [/?]  is  defined  as  the  partitioned  matrix 


au[B ] 

^12  \B\ 

■  B] 

®2l t-®] 

•  aln[B] 

fAB]  ■■ 

'  a,nn[B] 

(A.  1 ) 


[A]  ®  [i?]  is  of  order  ( mrxns ).  It  has  mn  blocks;  the  (/,  j)th  block  is  the  matrix  at-\B]  of  order 
(rxs). 
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B.  The  Mueller  Matrix  and  The  Kennaugh  Matrix 
The  Mueller  Matrix 

For  the  purely  coherent  case,  the  Mueller  matrix  [M]  can  formally  be  related  to  the  coherent  Jones  scattering 
matrix  [T]  as 


\m]=[ i  i  i  -1(m®[rr)[^ri=[^]([r]0[rr)[^r 

with  the  4x4  expansion  matrix  [A]  given  by: 
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so  that  the  elements  Mu  of  [M]  are: 
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If  [T]  is  normal,  i.e.  [T][T]r  =[7’]/  [T]  ,  then  [M]  is  also  normal,  i.e.  [MJfM]^  =[M]y[M] 


~iT* 


The  Kennaugh  Matrix 

Similarly,  for  the  purely  coherent  case,  [if]  can  formally  be  related  to  the  coherent  Sinclair  matrix  [.S']  with 
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[Ari=-[AT  as 


[if]  =  2[^f-1([5]0[5D[^r 


(B.3) 
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If  [5]  is  symmetric,  S  =  S  ,  then  [AT]  is  symmetric,  Ku  -  K ..  ,  so  that  for  the  symmetric  case 
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